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1.  Introduction 


The  objective  of  this  report  is  to  document  recent  one-dimensional  (1-D)  calculations  of 
temperature,  conductive  heat  flux,  and  heat  wave  velocities  for  an  opaque,  semi-infinite  cylinder 
due  to  radiant  heating.  The  Fourier  conductive  heat  flux  law  is  assumed  to  be  valid.  An 
idealized  model  is  used  for  the  calculations.  Radiant  flux  enters  the  planar  surface.  It  is 
uniform  and  steady  with  beam  and  cylinder  diameters  equal.  The  initial  temperature  is 
unifonn,  and  thennal  properties  are  constant.  Heat  wave  trajectories  (and  velocities)  are 
defined  by  the  value  of  a  constant  temperature  or  constant  conductive  heat  flux  that  propagates 
through  the  opaque  material.  For  the  constant  values  of  zero,  the  analytic  solutions  (1 )  of  the 
1-D  Fourier  conductive  heat  flux  law  and  transient  heat  transfer  equation  (derived  by  applying 
energy  conservation  to  the  Fourier  heat  flux  law)  lead  to  infinite  heat  wave  velocities.  For  the 
constant  values  >0,  the  numerical  solutions  at  ultra-short  transient  times  result  in  other 
nonphysical  velocities  (i.e.,  >3E(10)  cm/s). 

These  calculations  suggest  that  there  are  limitations  to  the  validity  of  the  Fourier  heat  flux  law 
and  heat  transfer  equation,  although  they  have  been  used  successfully  to  solve  engineering  and 
scientific  heat  conduction  problems  for  over  180  years. 


2.  Nomenclature 


Table  1  shows  the  units  and  values  of  the  thermophysical  properties  and  F  that  were  used  in  all 
calculations.  The  thermophysical  values  chosen  were  those  for  RDX  at  T  ~  288  K  (2). 


Table  1.  Nomenclature. 


Symbols 

Quantity 

Units 

Values 

t 

radiant  flux  time 

s 

— 

d 

distance  from  the  planar  surface 

cm 

— 

V 

temperature  increment,  i.e.,  V  (d,0)  =  0 

K 

— 

f 

conductive  heat  flux 

cal/cnT-s 

— 

X 

thermal  conductivity 

cal/cm-s-K 

2.5  x  10-4 

p 

density 

g/cm3 

1.76 

C 

heat  capacity 

cal/g-K 

0.3 

a 

thermal  diffrisivity  =  X/p-c 

cm7s 

4.7  x  10" 

F 

radiant  flux 

cal/cnf-s 

100 

1 


3.  Mathematical  Formulation 


For  the  simplified  model  with  opaque  (surface  absorption)  materials,  the  Fourier  heat  flux  law 


f(d,t)  =  -(l/X)-— V(d,t) 
Sd 


with  boundary  condition 


f(d,t)  =  F;  d  =  0,  t  >  0  . 


The  Fourier  heat  transfer  equation  is 


with  boundary  condition 


^V(d,t)  =  a|fv(d,t) 

dt  dd~ 


—  V(d,t)=-F/A,  ;  d  =  0,  t  >  0  . 
0d 


Using  equations  1  and  3,  Carslaw  and  Jaeger  (7)  show  that  the  Fourier  flux  transfer  equation  is 


with  the  boundary  condition 


~r  f(d,t)  =  a  f(d,t)  ;  d  >  0  ,  t  >  0 
dt  dd 


f(d,t)  =  F  ;  d  =  0 . 


The  solution  for  f(d,t)  (7)  is 


f(d,t)  =  F  •  erfc 


2- Va-t  J 


The  solution  for  V(d,t),  obtained  by  integration  of  equation  1,  is  (7) 

F  f  a-t  f-d2^|  d  f  d 

V (d,t) =  2 •  —  -  exp  —  ---erfc  — -j= 

X(  \  n  1  4at )  2  \2-Ja- 


Figures  1  and  2  are  calculated  values  of  f  (cal/cm2-s)  and  V  (K)  vs.  d  (cm)  at  t(s)  =  0.001,  0.01, 
0.05,  and  0.1. 
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Figure  2.  Temperature  vs.  distance. 

Average  heat  wave  velocities  for  conductive  flux  at  f  =  10  and  60  and  for  temperature  at  Y  =  50 
and  150  can  be  obtained  from  the  distances  between  intersections  of  their  corresponding  lines 
with  those  calculated  for  different  values  of  t.  Distances  vary  inversely  with  the  values  of  both 
f(d)  and  V(d),  which  suggests  that  the  largest  heat  wave  velocities  are  for  f  =  0  and  V  =  0.  To 
obtain  heat  wave  velocities  at  points  on  trajectories,  let  x  =  trajectory  distance  (cm)  from  the 
planar  surface. 

Values  of  x  at  given  values  of  t  can  be  obtained  from  numerical  solutions  of  equations  7  and  8. 
For  f  =  0  and  V  =  0,  the  numerical  solutions  (5)  are  unstable.  To  calculate  these  trajectories, 
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let  z  = 


2- y/a-t 


and  substitute  into  equations  7  and  8. 


V(z,  t)  =  2  •  — 


f  (z)  =  F  ■  erfc(z) . 

■{l-yfn  -z-ez  -erfc(z)). 


a-t 


1  1  it 

The  trajectories  (x  =  2  •  z  ■  yj  a-t)  can  be  calculated  for  given  values  of  t  from  solutions  of 
equations  7a  and  8  a  for  z. 

When  f(z)  =  0,  then  erfc(z)  =  0.  The  solution  to  equation  7a  is  z  =  go. 

When  V(z,t)  =  0,  then 

Vrtz  •  ez  •  erfc(z)  =  1 . 

Using  L’Hopital’s  rule,  when  z  =  oo,  then 

yfnz  ■  ez  ■  erfc(z)  =  1 . 

The  solution  to  equation  8a  is  z  =  go. 

For  V  >  0  and  f  >  0,  the  trajectories  (x  =  2  •  z  ■  yj a  ■  t)  can  be  calculated  from  the  values  of  x 
corresponding  to  given  values  of  t  from  numerical  solutions  of  equations  7  and  8  or  the 
corresponding  z  from  equations  7a  and  8a. 

Corresponding  velocities  —  x  (fv  for  flux  and  Vv  for  temperature)  can  be  calculated  by 

dt 

numerical  differentiation.  Calculating  trajectories  with  z  =  go  is  a  problem.  The  heat  wave 
velocities  can  also  be  obtained  analytically  from  the  calculated  trajectories  using  g(x,t)  = 


(7a) 

(8a) 


(9) 


(10) 


f 


constant :  —  x  =  - 
dt 


d 

ag(x’t}^ 


l 


-  g(x,t) 
dx  j 


When  f(x,t)  =  g(x,t),  then 


a 


fv  =  -^  =  z -J-. 
2  •  t 


(ID 


When  f(x,t)  =  0,  then  z  =  qo. 

When  x  =  0,  then  z  <  go  ,  i.e.,  0  <  x  <  go  . 

When  x  =  qo  and  t  =  go,  then  z  =  1  <  go,  i.e.  ,  0  <  t  <  go. 
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Using  equation  1 1  with  0  <  x  <  co  and  0  <  t  <  oo, 

fv  =  co. 

When  V(x,t)  =  g(x,t),  then 


(11a) 


4-t-a 


Vv  =  - 


n  ■  erfc 


x 

2-Vu 


n-z-ez  ■  erfc(z) 


(12) 


a  J 


When  V(x,t)  =  0,  then  z  =  oo.  Using  equation  10,  Vnz  •  ez  ■  erfc(z)  =  1 . 


Vv  =  z- 


When  0  <  x  <  oo  and  0  <  t  <  go,  then 


Vv  =  oo, 


(13) 


(13a) 


Equations  1  la  and  13a  show  that  when  0  <  x  <  go  and  0  <  t  <  oo,  and  f(z)  =  V(z,t)  =  0  ,  then 
Vv  =  fv  =  oo.  This  is  consistent  with  the  results  shown  in  figures  1  and  2. 


Figures  3  and  4  show  the  calculated  values  of  x(cm)  and  vV(cm/s)  vs.  t(s)  for  constant  values  of 
V.  The  x  ratios  (RxV)  and  vf  ratios  (RvV)  are  shown  on  the  right  y-axis.  In  figure  4,  the 
vertical  fine  indicates  that  for  V  =  1E(-19)  ;  vV  >  3E10  at  t  <1.37  10  23. 
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Figure  4.  Temperature  wave  velocity  vs.  time. 

Figures  5  and  6  correspond  to  figures  3  and  4  for  the  calculated  values  of  x  and  vf  vs.  time  for 
constant  values  of  f/F.  The  x  ratios  (Rxf)  and  vf  ratios  (Rvf)  are  shown  on  the  right  y-axis.  In 
figure  6,  the  vertical  lines  indicate  that  for  f7F  =  1E(— 10)  and  1E(— 1 6),  vf  >  3E(10)  at  t  < 

1.1  •  10  23  and  1.8  •  1 0  23,  respectively. 


Figure  5.  Flux  wave  trajectory  distance  vs.  time. 


1.284 


1.283 


1.282 


1.281 
1  28 


> 


- f/F  =  E(-16) 

- f/F  =  E(-10) 

- ratio 


Figure  6.  Flux  wave  velocity  vs.  time. 
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4.  Conclusion 


The  calculations  of  nonphysical  heat  wave  velocities  for  temperature  and  heat  flux  using  analytic 
solutions  of  the  Fourier  heat  flux  law  and  the  corresponding  transient  heat  transfer  equation 
indicate  that  there  are  limitations  to  their  validity.  This  problem  has  been  well  noted  in  the  heat 
transfer  literature  ( 4 ).  The  inability  of  the  Fourier  heat  transfer  (parabolic)  equation  to  predict 
measurements  of  heat  propagation  velocities  in  solids  near  zero  absolute  temperatures  and  for 
pulsed  laser  processing  of  solid  surfaces  involving  ultra-short  transient  heating  times  is  probably 
due  to  the  Fourier  flux  law  being  constitutive  and  not  valid  under  these  extreme  conditions.  The 
Fourier  equations  are  based  on  the  linear  dependence  of  heat  flux  and  temperature  gradient.  The 
details  of  the  mechanism  for  thermal  energy  transfer  (molecular  collisions  in  the  absence  of  net 
mass  motion)  are  still  under  investigation.  Consideration  of  phonon  interactions  as  a  molecular 
gas  (5)  has  recently  been  used  to  describe  the  thermal  energy  transfer  mechanism  for  conductive 
heat  flux. 

Non-Fourier-type  conductive  heat  flux  laws  in  solids  have  been  proposed  and  include  both  time 
and  temperature  gradient  dependence  (4).  These  laws  combined  with  energy  conservation  can 
lead  to  hyperbolic  type  transient  heat  transfer  equations,  which  avoid  the  Fourier  problem  of 
nonphysical  heat  propagation  velocity  predictions  and  give  better  agreement  with  measurements 
under  extreme  conditions.  Under  normal  conditions,  these  time-dependent  conductive  heat  flux 
laws  transit  into  the  Fourier  temperature  gradient  flux  law. 

Temperature  measurements  (using  different  techniques)  of  the  rate  at  which  equilibrium  is 
established  when  solids  (meat)  at  different  (near  room)  temperatures  are  placed  in  contact  have 
been  used  to  compare  heat  transfer  velocities  with  predictions  of  the  Fourier  and  non-Fourier 
type  conductive  heat  transfer  equations  in  solids  (6).  The  conclusions  as  to  which  type  equation 
is  in  better  agreement  with  the  data  are  in  disagreement.  It  is  of  interest  to  determine  the  lower 
time  limits  at  which  Fourier’s  heat  transfer  equation  remains  valid.  It  may  be  possible  to  obtain 
estimates  from  the  published  results  of  recent  pulsed  radiant  heating  experiments  on  opaque  and 
semi-transparent  solids. 
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